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We consider the discrete spectrum of the Dirichlet Laplacian on a manifold consisting 
of two adjacent parallel straight strips or planar layers coupled by a finite number 
iV of windows in the common boundary. If the windows are small enough, there 
is just one isolated eigenvalue. We find upper and lower asymptotic bounds on the 
gap between the eigenvalue and the essential spectrum in the planar well as 

for iV = 1 in three dimensions. Based on these results, we formulate a conjecture 
on the weak-coupling asymptotic behaviour of such bound states. 



1 Introduction 

There has been some interest recently to Laplacians on strips or layers. Such a 
system is trivial when the manifold is straight and the boundary conditions are 
translationally invariant, so there is a natural separation of variables. On the other 
hand, the spectral properties become nontrivial if the transverse modes are coupled, 
which can be achieved, e.g., if the manifold is bent, locally deformed, or coupled to 
another one |ES|, [DE|, PCTiSl gjTV| , [EVT|, |EV| . 
The interest stems from two sources. 



On the physical side, such operators with 
Dirichlet boundary conditions are used as models of various mesoscopic semicon- 
ductor structures. The corresponding solid-state literature is rather rich — see 
DEI, |ESTV|| for some references. On the other hand, bound states in systems with 



open geometries pose also mathematical questions such as the weak-coupling limit, 
validity of the semiclassical approximation, resonance scattering in such structures, 
etc. Some properties of them can be seen numerically [ESTV|| while analytical proofs 
are missing. Recall also that a closely related problem concerns Neumann Lapla- 
cians, namely the existence of trapped modes in acoustic waveguides [ELV| , pE|| . 
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In a recent paper [|EV1|| we have studied a pair of parallel Dirichlet strips of 
widths d\,d 2 coupled laterally through a window of a width 2a in the common 
boundary; we have shown that there are positive c\, c 2 such that the gap between 
the ground state and the threshold of the continuous spectrum can be estimated as 



c\a < e(a) 



< 



-c 2 a 



for any a small enough. The numerical result of |E^TV]| suggests that the true 
asymptotics is of the same type, but proving this assertion and finding the coefficient 
in the leading term remains an open problem. 

The aim of the present paper is to generalize the above inequalities to the case of 
a finite number of connecting windows and to a higher dimension. In the following 
section we shall prove the bounds for a pair of strips with N windows. In Section 3 
we formulate the analogous problem for two layers and prove two-sided asymptotic 
bounds for a single window shrinking to a point. Proofs rely in both cases on 
variational estimates and follow the same basic strategy as in [|EV1 ]. On the other 
hand, the existence of multiple windows or the change in dimension require numerous 
modifications, which prompts us to present the argument with enough details. 

The upper and lower asymptotics bounds we are going to derive are in each case 
of the same type differing just by values of the constants. We are convinced that 
ground state has an asymptotic expansion and its lowest-order is given by functions 
analogous to our bounds. This conjecture is formulated in the concluding section. 
At the same time, our present method does not allow to squeeze the bounds, or even 



to come close to the true values as the Remark 2.2 below illustrates. 



2 N windows in dimension two 

Consider a straight planar strip £ := M x [—d 2 ,di] • Given finite sequences C = 
{xk}k=i of mutually distinct points of the x-axis and A = {ak}k=i with > , we 
denote W4 := [xk— ak, Xk+^k] and set W := \Jk=i Wk ■ Then we define H(d\, d 2 ; W) 
as the Laplacian on £ 2 (£) subject to the Dirichlet condition at y = —d 2 , d\ as well 
as at the M \ W part of the this operator coincides with the Dirichlet 

Laplacian at the strip with the appropriate piecewise cut (see Fig. 1) defined in the 
standard way [|RS4 Sec. XIII. 15]. Following the notation introduced in UEVlH we 



put d := max{c?i,(i 2 } and D := di+d 2 . If d\ = d 2 , the operator decomposes 
into an orthogonal sum with respect to the y-parity; the nontrivial part is unitarily 
equivalent to the Laplacian on L 2 (S + ) , where S + := M x [0, d] , with the Neumann 
condition at window part W of the x-axis and Dirichlet at the remaining part 
of the boundary; we denote it by H(d; W) . If the specification is clear from the 
context, we will often denote the operator in question simply as H . 

We need a quantity to express the "smallness" of the window set. We define 

N N 

J(W) := Y, a k\m\ = 2$> 2 ; (2.2) 

k=l k=l 
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Figure 1: Window-coupled planar waveguides 



then the result of |Wj]] generalizes to the present situation as follows: 

Theorem 2.1 a ess (H(di, d 2) W)) = [(n/d) 2 oo). The discrete spectrum is contained 
in ((n/D) 2 (rr/d) 2 ) , finite, and nonempty provided W^l // ^ s sufficiently 

small, o~$ lsc {H{di, d 2 \ W)) consists of just one simple eigenvalue e(W) < {^/d) 2 and 
there are positive ci, c 2 such that 



-rJ(W) 2 < e(a)-(^) < -c 2 I(W) 2 



(2-3) 



holds for any /(W) small enough. 



Proof: (a) The upper bound. In the symmetric case, d\ = d 2 , the trial function will 
be chosen as ip = F + G , where 



with 
and 

with 



F{x,y) := fi{x)xi{y) , 
h{x) := m a x{x [xi - ai , XN+ a N] (x), e^"^, e -«l-«*-*l} 



(2.4) 



JV 



:= ^G k {x,y) 



k=l 



\W\ 



X[x 



k —a k ,x k +a k 



(*) cos (^-^) 



(2.5) 
(2.6) 



where | W| := 2J2k=i a k , and 

e -iry/2a k 



R k {y) ■■- 



2(1- I) e-^/ 4 ^ 



1/ e 



o - 



2' 



(2.7) 



for fc = 1, 2, . . . , N . As before Xn{y) = \f\ sin [^f] , n = 1, 2, . . . , denote the 
"transverse" eigenfunctions — to be not confused with the indicator function xm of 
a set M . Notice that as long as we work with trial functions of Q{H) , the window 



smoothing employed in |EV1| is in fact not needed — cf. 
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The functional L(ip) := (Hip,ip) — (X\ H^ll 2 can be expressed as 

/ 7r \2 / 9 N „ x . fc _|_ afc 

= ll^f+II^H 2 - - \\G\\ 2 -2-J-Y: / G k (x,0)dx. (2.8) 
\a/ a \ a ,_, j Xk -a k 



Since f x , G x have disjoint supports, we have H^H 2 = ||F a ,|| 2 + J2k=i \\Gk,x\\ 2 , where 
Gk,x '■= d x Gk ■ The fcth term of the last sum equals rftjx 2 ay\ W|~ 2 ||-Rfc|| 2 2( d ) , and 

\\Mh m = -+(i-^) e^<^(l + e 1 ) 

y ' ix yo IX J IX 

for any E\ > and a k small enough. Obviously, Q-au Gk( x , 0) dx — ^ r^afj Wl -1 , 
and furthermore, a bound to ||Gfc, y || 2 follows from 

\\B' II 2 , - — + f- - J!^ P -^/2« fc < 

II^IIlw) - 4afc + ^ 4a J e < 4afc 

for a fc < 7rd/8 , which means that ||(j/fc l2/ || 2 < tt J2k Vk a t\ W|~ 2 . Now we can put these 
estimates together using ||-F X .|| 2 = k ; neglecting the negative term — ( UGH 2 , we 
arrive at the inequality 

t(,u\ ^ * 16y/2 ^ Vk a 2 k * ry 2 q 2 

iW<K -^£w + 7r(2+£l) Sw 

The sum of the last two terms at the r/is is minimized by — 7rd 3(2+ ei ) a fc • To 
conclude the argument, we have to estimate the trial function norm H^H 2 from 
below. The tail part is HVil^eMW = while the window contributes by 

N 2 3 

IH| 2 6 w < 2 \\ F \\lew + 2 \\ G \\leW = \x N -x x +a N +a x \ + 4 £ \\R k \\l2 m , 

so H^H 2 > (1 — ^2)^ 1 holds for any e 2 > provided |W| is small enough. Mini- 
mizing the obtained estimate of L( / 0)/|| , 0|| 2 over we find 

< - l- £2 -i (2.9) 
W 3 (2 + £l ) ^ / 



which yields the upper bound in (|2.3| ) for d\ = di . The extension to the nonsym- 
metric case proceeds as for N — 1 ; the trial function is chosen in the above form 
for the wider channel, while in the narrower one it is given by (|2.5| ) transversally 
rescaled. 

Remark 2.2 The bound can be improved, for instance, by replacing the factorized 
form (|2.6| ) by a series, whose terms will be products of the trigonometric basis in the 
window with the functions Rk, n {y) decaying as exp{ — around y = (in the 
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above estimate we used just the first term of such a series). However, the gain is not 
large. To illustrate this fact, take N = 1 and d = it . The use of the series leads 
then to the upper bound (^j improving the coefficient by (7r 2 /8) 2 ~ 1.52. A 

comparison to the numerically determined ground state fESTVf shows that the true 
asymptotic behaviour should be ~ (2.23 a) 4 , so the obtained C2 is still two orders 
of magnitude off mark. The reason is obviously that the wavefunction is affected by 
the window outside the transverse "window strip" as well. 

Before proceeding to the lower bound, let us state some auxiliary results: 

Lemma 2.3 Let J[0] := f b a {<f>'{t) 2 + m 2 0(t) 2 ) dt for E C 2 (a,b) with 0(a) = c a 
(a fixed number). Given niQ > ; there is «o > such that 

J[(j)} > a mc 2 a (2.10) 

holds for all m > mo . 

Proof: The mimimum is obviously reached with (j)'(b) = . The corresponding Eu- 
ler's equation is solved by 4>o(t) = die~ mt +d2e mt , where d\ = c a {e~ ma + e m ( a ~ 2b )^ 
and c?2 = die~ 2mb . Since m _1 c~ 2 inf J((f>) > for any m > mo, it is sufficient to 
check that fl2T0D remains valid as m — > oo ; evaluating the functional for 0o we 
find lim^oo J (4>) = mc 2 a . | 

Lemma 2.4 Suppose that minimizes J[0] : = J 2a (<f>'(t) 2 + p 2 (t>{t) 2 ) dt for positive 
a, p within C 2 (a,2a) with the boundary condition 0(a) = c a ; then 

|0(2a)| < 2\c a \e- pa . (2.11) 

Proof: Assume for definiteness that c a > . By the mentioned symmetry argument 
again, 0'(2a) = 0, and its explicit form is <f>(t) = c a coshp(2a — t)/ cosh pa, which 
yields 0(2a) < 2c a e~ pa . I 

For the sake of completeness we reproduce also the following assertion the proof 
of which is given in pVl]| : 

Lemma 2.5 Let e C 2 [0,d] with 0(0) = (3 and 0(d) = . If ((f), Xi) = 0, then 
for every m > there is do > such that 

(%'{t) 2 dt + f-V T0(t) 2 dt - f^V (Uitfdt > ^! (2.12) 
Jo \ a / Jo V a / Jo a 

holds for all a small enough. 



(b) Proof of Theorem \2.\ continued: The lower bound is again the more difficult; 
however, we may restrict ourselves to the symmetric case only because inserting an 
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additional Neumann boundary into the window we get a lower bound, and therefore 
we consider in the following the spectrum of H = H(d; W) . 

We begin with a simple observation that it is sufficient to estimate L(ip) := 

(Hip,ifj) — (jj) \\ip\\ 2 from below for all real ip of a core of H, say, all C 2 -smooth 
ip G L 2 (S + ) satisfying the boundary conditions, since H commutes with complex 
conjugation. The main difficulty brought by the existence of multiple windows is 
that we are no longer allowed to restrict ourselves to trial functions symmetric with 
respect to the window centers. The strategy we employ is to split from the beginning 
a part of the kinetic-energy contribution to the functional, say, | HV^II 2 ? which will 
be at the end used to mend the problems coming from the asymmetry, i.e., we begin 
with estimating Lq{jP) := L{ip) — | H^H 2 . 

A trial function of the indicated set will be written in the form of a Fourier 
series, 

oo 

1>(*,V) = J2 C n( X )Xn(y) (2.13) 
n=l 

with smooth coefficients c n (x) = (ip(x, ■), Xn) , which is uniformly convergent outside 
the windows, x £ W . We split further the lowest transverse-mode coefficient by 
putting 

N 

/i:=ci-£A> (2.14) 



k=l 



where 



p | c k (x)-a k ... x G [x k -2a k ,x k +a k ] (2 15) 

... otherwise ^ ' ' 



with a k := c\{x k — 2a k ) , i.e., each one of the functions f k vanishes at the left 
endpoint of the appropriate extended window; in contrast to | EV1 | we double the 
left half of the window only. Writing the full trial function as 

V(x, y) = F(x, y) + G(x, y) , F(x, y) := fi{x) X i(y) , (2.16) 

we can cast the reduced energy functional into the form 

3 /7T\ 2 N 7T F> r 

MV) = i ll^f + ||Gj 2 - (-) IIG|| 2 -E 2ak duL k G(x ' 0)dx - (2 ' 17) 

Contributions to (|2.17 ) from different parts of the strip S + will be estimated sep- 
arately. The out-of-window part consists of the sets 



£*>i = {{x,y) 
u k = {(x,y) 
u N+l = {(x,y) 

The expansion ( |2.13[) yields 



x < Xi — ai } , 

Xk-i + a^x < x < x k -a k } , k = 2,...,N. 
x > x N + a N } . 
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1 00 r °°/-7r\2 r 

4 „ =1 -^fc „=iV«/ v J ^fe 



and therefore 



ML fe + II^IL fc 



with some fiQ > follows from Lemma [O] (applied to c n (—x) ) for k = 2, . . . , N . 
The same inequality for fc = 1 is derived as in |[EV1|| ; for the right tail we use just 
the fact that the expression is positive so we can neglect it. Since ip x = G x inside 
the (left extended) windows, we arrive at the bound 



N 

E 

k=l 



-u \\r ii 2 



_ iir ii 2 

4 W^XUX^Hi, 



- 'd) 



\\G\\xeW k 



(2.18) 



00 7T [2 f ~] 

+ no J2nc n (x k -a k ) 2 - 2a k -\ - / 0) dx \ 

^2 d \ d Jw k J 



Our next goal is to estimate the contribution to ||Gz|| 2 from the extended windows, 
&k '■= [%k— 2afc, Xk+dk] ■ In distinction to the case N = 1 , however, even the lowest- 
mode projection of G may not vanish at the right endpoints of these intervals, 
so the inequality (5.6) of [ EVI | has to be modified. Fortunately, it is sufficient to 
change the coefficient: if a function G : S + — ► C 2 (S + ) vanishes for x = x k — 2a k , 
the inequality (4.2) of [EV1] in combination with a symmetry argument imply 



1 a, 



x\\xee k 



> 



(-)' 



\r\\ 2 



(2.19) 



To use this result we split the function by singling out the projection of G onto the 
first transverse mode, 



N 



G(x,y) = G 1 (x,y) + G 2 {x,y), G 1 {x,y) = £ f k {x) X i(y) ■ (2.20) 

k=l 



We have 
1 „ 



1 



1 



1 



\x<=£ k \W k 



x\\xeWk 



and therefore 

1 N f 1 /7T\ 2 

LoW > 2 \M\l?s + £ {jH^fc + H^vll^w* " {2) W G Wlew h 

[2 



(2.21) 



71 2 f 

+ /i V nc n (x k -a k ) 2 - 2a k - \ - / G(x, 0) dx + 
n=2 a V « -'Wfc 



1 /jr_y 

2 \6aJ 



G i\\le£ k } 
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with £ : = UfcLi^fc- To proceed further we split the function G2 in the k-th 
extended window as G2(x,y) = G(x,y) + F(x,y) , where 

00 

r(z,y) := c n (x k -2a k )xn(y) ■ 

n=2 

The second part is independent of x while the first one vanishes at left endpoint, 
so G 2 , x = G x may be estimated by means of (|2.19|) and the Schwarz inequality as 

2 



k 



(2.22) 

2 _( M 2 11^112 
2 V6a fc 7 " *" n * V6a fc / 



where we have denoted Q k '■= &k x [0, a&] • To make use of the last estimate we have 
to find an upper bound to ||r||^ . To this end we notice that 



r2 



(i) instead of assuming c n G C , the lower bound can be looked for in a wider class 

of ip with piecewise continuous coefficients, 

(ii) on the other hand, we may restrict ourselves to those ip which satisfy for 
x G £k \ Wfc and n>2 the inequality 



cosh (5 \/n 2 — 1 (x — Xk + 2a,k 

\c n (x)\ < C(x) := \c n (a)\ W , ^ . (2.23) 

cosh ( v n — 1 



To see that we split the trial function in analogy with [EVT 
#(x,y) : = 



ip(x, y) - c n (x)xn(y) ■■■ x eS k \yVk 
i[>(x, y) ... otherwise 



The basic expression L(?/>)/||?/>|| 2 can be then rewritten as 



dx 



+ J2k=l Is k \Wk C n( x ) 2 dx , 

where L(V>) := / s+ (J^l 2 + iV-vl 2 ) (x,y) dx dy . We may assume only those 
■0 for which the numerator is negative; the part of its last term correspond- 
ing to the "window neighborhoods" is minimized by the hyperbolic function 
c^f of ( |2.23j ) (see the proof of Lemma |2.4| ). It follows that replacing c n (x) 2 



by mm{c n (x) 2 , c*^ (x) 2 } we can only get a larger negative number, while the 
positive denominator can be only diminished. 



8 



To estimate the norm of T restricted to Q k , we adapt again the argument of |[EV1| 
and divide the series into parts referring to small and large values of y , and employ, 
respectively, the smallness of || Xtx |^ [0, a.] || and the bound ( [2.23| ). This yields 

/no \ 2 



^ dx J q dy [ c n [2a k }xn(y) 



\n=2 



^2 c n [2a k ]xn(y) | dy + Qa k 

V 



n=2 



c n [2a k }xn{y) I dy 

\2<n=[a~ 1 }+2 



< 24a fc | n ' C «N 2 / Xn(y) 2 dy ^ n 

n=2 J ° J y n=2 



I 



+ 24a k 



nc n [a k f J Xn{yfdy\\ 

\2<n=[a~ x ]+2 ' " / \2<n=[a- 1 ]+2 



n -l e -(27ra fc /d) N A?3T 



where c n [ja k ] := c n (x k —ja k ) and [■] denotes the entire part; in the the last step we 
have used the bound |c n [2a/-]| < 2|c n [a^]| exp |— ^y/n 2 — l| which follows from 
Lemma 2.4 . In analogy with [EV1], this implies the existence of a positive C k such 
that 

oo 

\2 



|r||n fc < C k a\ J2 nc n( x k-ak) 



(2.24) 



n=2 



From now on we consider continuous coefficient functions again. By ( p. 22 ) we have 
1 



\ G 2,x\\lte k + Y nc n( x k-akf 



> 5 



71 



I2a u 



n=2 



l G 2||L 



5 ( 7T 



2 \6a k J 



\ 2 oo 

) \\ T \\n k + Mo Y nc ni x k-ak) 



n=2 



for an arbitrary 5 G (0, 1] ; if we choose the latter sufficiently small, the sum of the 
last two terms is nonnegative for each k = 1, . . . , N due to ( p.24| ), so 



N 



+ £ \\\g 



"Turn 2 
d) 11 lUw « 



2 

2|in fc 



(2.25) 



1 / n 

2 \6a k , 



xes k ( > 



7T /2 /" 

— 2a k — \ - / G(x, 0) dx 
a V « w fc 

where we have denoted m := £ \/#. 

Next we express the first term in the curly bracket using the decomposition 
(|2.20D , properties of the transverse base, and an integration by parts, 



\r t II 2 _ wn n2 



\G 



2,y\\xeW k 



7T 2 

d \ d 



f k (x)G(x, 0) dx. 



As in |[EV1|| we estimate the last term by the Schwarz inequality, substitute into 
( |2.25| ), neglect ||Gi,i,||x e yv fc as well as 

^ u r ,,2 7T(7T + V^) 2 



72a 2 ll-HI^ rf 2 W^VWXEW, 

which is positive for a k small enough, obtaining 



1 N f 

z k=i 1 



(2.26) 



^l|G(-,0)|| 2 eWfc -2«^w| / G(*,0)dx 



)V, 



By Lemma |2.5|, the sum of the first three terms in the curly bracket is bounded from 
below by i\\G(;0)\\l 

ew k f° r some d k > . Since (d k /2a k ) — (irV2/d) > holds 
for a k small enough, we have 



LoW > ^11^11^+ E{^l|G 2 (-,0)|| 2 eWfc -4a^^||G 2 (-,0)|U. eWfc }, 



where we have employed again the Schwarz inequality. The k-th term of the sum 
reaches its minimum w.r.t. the norm at — J^3 af-Ci 2 ,. Returning to the original 
functional and neglecting in the first term of the last estimate all contributions 
except the one coming from the leftmost component of M \ £ , we see that there is 
a positive 7 such that 

1 1 N 

U$) > 4 + g W^Wl<* k -2a k - 7 E «M (2-27) 



k=l 



holds provided |W| is small enough. 

To conclude the proof, we denote £ k : = x k —2a k —x\+2ai and employ the identity 

N N N 

E <*A = E <*W k + E(«*- a >l ( 2 - 28 ) 
*!=i fe=i fe=i 

together with the estimate 

-\M\ 2 > -KW 2 > — E { - k ~ ai)2 

4 ll</yi11 - 4 11 111 - 4AT ^ 

If -y(al- aj)a 2 k + ^(a k - «i) 2 > holds for all k = 2 . . . , N , the bound ( ^27|) 
reduces to 

1 * 

£(V0 > 2 ll^H^r**- ^1 E a *- ( 2 - 29 ) 

fc=l 
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On the other hand, suppose that the endpoint values satisfy — a.\ = 0(ak) as 
a k -> for k E JC (Z {2, . . . , N} . In view of ( ggjp we have 



* f 1 1 

- k=i fce/c ^ 



however, the last term is O (J2k&tc a t) j so ( |2.29| ) is valid again with a smaller positive 
coefficient in the last term. Since ||^|| 2 > 2 /f^ 2 " 1 c\{x) 2 dx , the quantity of interest 
is bounded from below by 

> - rxi-2a 



2/« Cl (x) 2 dx 

The rhs is minimized by the function c\(x) = aie K ^ x ~ Xl+2ai ^ which yields the value 
(k 2 /2) — 7_T(>V)k ; taking the minimum over k we find 

^| > - 7 2/(W) 2 . | (2.30) 



3 Window— coupled layers 

The setting of the three-dimensional problem is similar. We have a straight layer, 
S := M 2 x [— d2,d\] , and a set W C M 2 which can be written as a finite union, 
W := U^ 1 Wfc , whose components are open, connected sets of nonzero Lebesgue 
measure; without loss of generality we may suppose they are mutually disjoint. Then 
we define H(dx, d 2 ; W) as the Laplacian on L 2 (Y,) obeying the Dirichlet condition 
at the boundary of E, i.e., y — —d 2 ,di, as well as at M 2 \ W. This operator 
coincides again with the Dirichlet Laplacian ||RS4|, Sec. XIII. 15] for the sliced layer 



the two parts of which are connected through the window set W . We use the same 
notation as above, d := max{di, d 2 } and D := d\ + d 2 . The nontrivial part of the 
symmetric case, d\ = d 2 , reduces again to analysis of the Laplacian L 2 (E + ) , where 
E + := M 2 x [0, d] , with the Neumann condition at window part of the plane y = 
and Dirichlet at the remaining part of the boundary; this operator will be denoted 
as by H(d; W) . 

Our main aim here is to prove a weak-coupling asymptotic estimate for a pair 
of layers connected by a single window. 

Theorem 3.1 a ess (H(di, d 2 ; W)) = [(ir/d) 2 , oo). The discrete spectrum is contained 
in ((tt/D) 2 , (n/d) 2 ) , finite, and nonempty provided W / I. Suppose further that 
N = 1 and W = aM for an nonempty open set M contained in the unit ball 
Bi C M 2 . Then crd isc (H(di, d 2 \ aM)) consists of just one simple eigenvalue e(aM) < 
(7t/g?) 2 for all a small enough, and there are positive ci, c 2 such that 

- exp (-cigT 3 ) < e(a) - (^j < - exp (-c 2 cT 3 ) . (3.1) 
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Proof is based again on variational estimates. The upper bound in the symmetric 
case, di = d 2 , employs the trial function ip = F + r]G , where F(x, y) := fi(x)xi(y) 
again with 

fi(x) := mmi 1, — — — -} , (3.2) 



and 



G(x,y) ■= XaM(x)M x ) R (y) 



(3.3) 

where 0^ is the ground-state eigenfunction, \[<j>\ || = 1 , of the operator — A" M 
corresponding to the positive eigenvalue /ii(a) = ^i(l)a -2 , and 

e -\f^l{ a )y _ _ _ y £ 



H(y) 



- 



(3.4) 



2(l-jf)«p (-ivW 



2/ e 



2' 



d 



Using — Xi = (7r/d) 2 xi , a simple integration by parts, and the fact that the vector 
functions V/i and have disjoint supports, we can express the reduced energy 

as 



functional L(if)) := (i*V,V0 - (§) 

£(V0 = ||V/i||i 2(JR2) + r/ 2 ( /ix(a) - (~ 



L 2 (0,d) 



v 2 \\R'\\h m - 2 vx'M I __<Pr(x)dx 



ahl 



(3.5) 



where the negative term in the bracket can be, of course, neglected. The second and 



the third term at the rhs can be estimated in analogy with |EV1 



^i(a)||-R||z2 (0)d) -77 || B!\\ L 2 { q 4) < 



2a 



(2+ei) 



for a fixed £1 > and any a small enough. In a similar way, the last term equals 
— 2t]Xi(0)Ca , where C := j* M </>i (a;) dx . Finally, the first one can be evaluated by 
means of |A|, 9.6.26], gBM] , 1.12.3.2], 



^oM 2 ||VA||l 2 



(jp) 



2tt 



i /t 2 a 2 K;(/ta) 2 - - (kV+ l) ^(ko) 



Using — K[(£) = K (£) + £ 1 Ki(£) in combination with the asymptotic expressions 
Ko(£) = -ln£ + 0(l), ^1(0=^+0(1110, we find 



II Williams) < 



2tt(1 + £ 2 ) 
ln/ca 



for a fixed e 2 and a small enough. Substituting these estimates into ( |3.5| ) and 
taking a minimum over 77 we arrive at the bound 



2tt(1 + £ 2 ) 
In Ka 



2x'i(0) 2 C 2 



(2 + £l W Ml (l 



(3.6) 
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It remains to find a lower bound to 



> 



\x\>a ~ ^\\ r \\\x\<a ~ 2?? 2 ||-F||f, c |<a — UV^ll |a;|>a 



-2IIFI 



- 2na 2 



2r] 2 \\R u2 



L 2 (0,d) ■ 



The last term is 0(a) , while the first one can be expressed as 

7T 



K 0{^o) 2 \\F\\ 2 \x\>a 



Ki(ko) — Kq(ko) 



K 



+ 0(a 2 \nKa) ; 



using the asymptotic behaviour of K we find 
fixed £3 > and a small enough. Hence 



> txk 2 (ln/ta) 2 (1 — £3) for a 



< 



k 2 In KCL 
vr(l-e 3 ) 



(Da 3 In na + E), 



(3.7) 



where := 27r(l + e 2 ) and 



2 X [(0) 2 C 2 



(2 + 50^(1) 



Minimizing the r/is of ( |3.7| ) with respect to k , we conclude that to fixed positive 
£1, e 2 and £3 G (0, 1) there is a function g such that 



L(V>) 



and 



0(a) 



1 + g 2 J_ e -2E/Da 3 

I — £3 a 2 



(3.8) 



as a — ► . The upper bound in ( |3.1|) follows readily from ( p.8|) ; the extension to the 



nonsymmetric case is obtained as in [EV1 



Remark 3.2 In fact, one could suppose M = B\ because the eigenvalue is pushed 
up if we reduce the window to a circle contained in M , and the obtained bound 
is all the same not optimal as in Remark [2^. In the rest of the proof we embedd 
M into a circle leaving the question about relations between the constants and the 
geometry of M to more sophisticated methods. 

The lower bound can again be proven in the symmetric case only. We begin with 
auxiliary results. When constructing the trial function component ( |3.2| ), we have 
used implicitly the fact that the functional F : F(4>) = J a °° (0'(t) 2 + m 2 <j)(t) 2 ) tdt on 
C 2 ([a, 00)) with the condition <p(a) = a and fixed positive a,m and is minimized 
by 

K (mt) 



0o : 4>o(t) = a 



Ko(ma) 



(3.9) 



as can be easily seen from solution of the appropriate Euler's equation. Furthermore, 
a two-dimensional analogy of the bound (4.2) in [EV1| is given by the Friedrichs 
inequality | [Ne| , Thm. 1.9]: if Q C M n , n > 2 , is a bounded domain with Lipschitz 
boundary, there is a positive c such that 



iiv/ir> C 



(3.10) 
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holds for every / G (Q) . The constant is, of course, easy to find for the circle 
Q = B a in terms of the appropriate Bessel zero, c = j^a -2 . 

Repeating the argument of [|EV1|| and the previous section, we infer that one has 



to find a lower bound to L(ip) / \\if)^ 2 over all real if) G L 2 (E) , which are C 2 , radially 
symmetric, and vanish at the boundary except in the window. We can express such 
a if) in the form of the series ([2.13|) again, where the convergence is uniform for 
\x\ > a. The coefficients c n depend in fact only of r := \x\ . Moreover, in analogy 
with ( |2.23|) we may restrict our attention to trial functions with 



K U s/tf=l> 
\cn(r)\ < \cn(a)\ Y (3.11) 
K Q (J Vn 2 - la] 



for n > 2 . As before we introduce 

F(x,y) := { a f[ V \ , °^ 2a (3.12) 

y ' \ c x {r)x\{y) ... r>2a K J 

with a := ci(2a) , and divide the rest G(x,y) = ip(x,y) — F(x,y) into 

Gi(x,y) := (ci(r) - a)xi(y) 

supported in the extended window region, r < 2a , and G2(x, y) = G(x, y)+T(x, y) 
with 

oo 

Tfay) := c «( 2a )Xn(2/) • 

n=2 

We start estimating the reduced energy functional 

L(V>) = || VM 2 + \\Gy\\ 2 - Q)' \\G\\ 2 - 2^(0) jf G(x, 0) dx (3.13) 
from the "external" contribution to the first "two and a half terms, 



1 , ,,o ,, „ ,,o / 7T 



In ■= §l|V^I|j!>a+l|GX> a -K) ||G|| 2 > a 



2 



E jf (<(0 2 + 2 Q) 2 (n 2 - l)c n (r) 2 ) rrfr 



oo 

7T 

n=l 



> ^E/ a °(4W 2 +2(^) 2 Cn(r) 2 ) raV 



(3.14) 



> ii Vc n (a) —a ) f > — — Vnc„(a) 



where in the last line we have used ( |3.9| ), evaluated the integral as in the first part 
of the proof, and employed the inequality > ^o(0 which follows from the 

well-known integral representation [AS, 9.6.24]. Next we turn to 



VMU2a = IIV^iH^ + || V x G 2 || 2 <2a . (3.15) 
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By assumption, G\ vanishes at r = 2a , so the first term can be estimated from 
(ETiq) as 

II ^^l|lr<2a > ^ ll^l|lr<2a = ~jll^l|| 2 > (3.16) 

where 4Ci := jgi- Furthermore, introducing the window neighbourhood Q a := 
i?2a x [0, a] , we have 

II V a .G2||r<2a = II ^ xG || 2 < 2a > — ^ H^Hr^a 

(3.17) 

> SII^IIL > ^II^IIL > ^IIg 2 ■ 



Cl ||A||2 ^ ||/*t||2 ^ Hrr ||2 ||r||2 

a 2 IMIn 2o > fl2 ll^lln 2a > 2a2 \\^2\\n 2a fl2 INIn 2a 



for all a < d and 5 G (0,1]. The last norm can be estimated as in the previous cases 
by combining the smallness of the \ n norm restricted to [0, a] with the dominated 
decay (|3.11|) , 



|r||L = 4vra 2 / ]T c n (2a) X n(y) dy 

•' \n=2 / 

< Kvra 2 I " _1 Cn(a) 2 / Xn(y) 2 dy J2 n 

n=2 J ° J n=2 

( ra \ ^2^^231 

+ 87ra L " c>) X Xn(y)dl/ L ^.. ^F7^^ i 



v 2<n=[a- 1 ]+2 / 2<n: 



^.i ]+2 n^ 2 (f 



a) 2 



" d ^ + 2 ^_ i]+2 ^2 (fv ^ )y |^^ 
The sum in the bracket can be estimated as 

2<n=^- 1]+ 2 < (if V^I) ~ W (f V/F^T) " * (g) 

for a < |\/3 , and the integral on the r/is is convergent, because i^o(£) ~ \f& e ^ 



2? 

as £ — ► oo . Hence there is a positive C*2 independent of ip and a such that 

4ll r HL < C 2 aJ2nc n (a) 2 . (3.18) 

a n=2 

Combining the estimates ( p.l4| )- (|3.18| ), we arrive at 

( 2 \ oo /-» 

- 5C 2 j ^ncniaf +^\\Gi\\ 2 +^ \\G 2 \\l a , 
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which gives 



Cx 

2 



a 



Li+ L 2 > —\\Gi\\ 2 -\ \\G 2 1| q 



(3.19) 



for some m > and all sufficiently small a . 

in 1 

2tt 



The norm of G y is estimated as in the two-dimensional case |[EV 1 



|^y||r<a — ll^2,j/||r<a 



d 2 



(2||G 1 ||?< a + rf||G 2 (-,0)|| r 2 < a 



which together with ( |3.19| ) yields 



Lx + L 2 



\r ii 2 



\r\\ 2 



> \\G " 2 



2,J/llr<a 



+ ^rWG 



2lln« 



C 27T 



7T 



l^2||r<a 



|G 2 (-,o)|| 2 <a 



Co 



'<7l _ 7T(7T + 4) ' 

a 2 d 2 



lr> II 2 

1^1 llr<a 



> - - -T HG 2 (-,0)|| 2 <a > ^||G 2 (-,0)|| 2 <a 



2a 



for a positive cq and any a small enough; in the second step we have neglected a 
positive term and employed Lemma p~5| . Substituting from here to ( |3.13| ) and using 
the Schwarz inequality, 



/ G(x,0)dx < \\G 2 (;0)\\ r<a ^a 

JBn 



we get 



> ^liv.vil 



2aa X , 1 (0) v ^||G 2 (-,0)|| r < a + ^||G 2 (-,0)|| 2 <a 



r>2a 



2 

r>2a 



2a 



2™V 1 (0) 2 3 
c 



The first term on the rhs can be estimated from below by the first transverse-mode 
contribution. The same applies to ||^|| 2 , so finally we find 



m > ^ a c\{r) 2 rdr-^a* Cl {2af 
2 J£ d(r) 2 r dr 



(3.20) 



In analogy with ( |3.9| ) one has to solve the appropriate Euler's equation to check that 
the rhs of ( |3.20|) is minimized by C\ = <ft K for some k > , where 



(j> K (r) ■= cx{2a) 



Ko(2na) 
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Substituting into Q3.20| ), evaluating the integrals, and taking the asymptotics for 
small a , we infer that 

LW > -nHn(2na) ( n ^ ' q» \nj2n a) ' 1 ~ ^ 



holds for any fixed e±, e<i > and all sufficiently small a . It remains to find the 



minimum of the rhs with respect to k . However, since it differs from ( |3.7|) just by 
the values of the constants, the argument is concluded as in the first part of the 
proof. I 



4 Conclusions 

To make sense of the derived bounds one has to take into account two aspects of 
the problem. First of all, we have mentioned already that the discrete spectrum 
can also be found numerically by means of the mode-matching method; a detailed 



description of the two-dimensional case is given in |[ESTV| . Although the method 
converges rather slowly if the window is narrow, the results obtained for a single 
window clearly suggest that the true asymptotics exists and is of the same type as 
our asymptotic bounds. 

Another insight can be obtained from comparing our result with the well-known 
weak-coupling asymptotics for Schrodinger operators in dimension one and two 



BGS| , [Kj , . The ground state of the coupled strips in the narrow-window case is 



dominated the lowest transverse-mode component with long exponentially decay- 
ing tails and a local modification in the coupling region. In a similar way, a link 
can be made between window-connected layers and a two-dimensional Schrodinger 
operator. The comparison shows that the attractive interaction due to opening a 
narrow window (in particular, by changing the Dirichlet b.c. to Neumann at a short 
segment of the boundary in the symmetric case) acts effectively as a potential well 
of a depth proportional to the size of the window. 

Conjecture 4.1 Let H(dx,d2',W) be the operators described above. The ground- 
state eigenvalue behaves for small \W\ as 

e(a) » (£f - I (£c 2 ,^k) ... dimS = 2 (4.1) 



e (a) 



7T\ 2 1 f 



d) -exp{-[J2c 3 , k ^y k ) \ ••• dimS = 3 (4.2) 



\fc=i 



where v :— d 1 minjc?!, <i 2 } , and a k in the three-dimensional case is the scaling 
parameter of the k-th window. 

The conjecture is based on the described analogy only, and therefore it is difficult to 
say more about the coefficients. It is not excluded that they depend on the geometry 



17 



of the window-center set for N > 1 ; in the three-dimensional case the shapes of the 
scaled windows may also play role. We refrain from speculating about the nature of 
the error terms. 

On the other hand, we are convinced that the open "constant cross-section" 
shape of our regions E is crucial for the asymptotics. For instance, if E is instead 
a bounded planar region with the Dirichlet boundary in which we open a window 
(to another bounded region the essential spectrum threshold of which is not lower) 
or a Neumann segment, we conjecture that leading term in the ground state shift is 
proportional to the square of the window width. Moreover, the same asymptotics is 
expected to be valid for higher eigenvalues provided the corresponding eigenfunctions 
are locally symmetric with respect to the window axis. In any case, proving of such 
asymptotic properties represents an intriguing mathematical problem. 
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